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Abstract 

We argue that quantum Liouville field theory supplemented with a suitable source term is the effec- 
tive theory which describes the short-range correlations of the gluon saturation momentum in the two- 
dimensional impact-parameter space, at sufficiently high energy and for a large number of colors. This 
is motivated by recent developments concerning the stochastic aspects of the high-energy evolution in 
QCD, together with the manifest scale invariance of the respective evolution equations and general con- 
siderations on the uncertainty principle. The source term explicitly breaks down the conformal symmetry 
of the (pure) Liouville action, thus introducing a physical mass scale in the problem which is identified 
with the average saturation momentum. We construct this source term for the case of a homogeneous dis- 
tribution and show that this leads to an interesting theory: the relevant correlation functions are ultraviolet 
finite (and not just renormalizable) when computed in perturbation theory, due to mutual cancellations of 
the tadpole divergences. Possible generalizations to inhomogeneous source terms are briefly discussed. 



E-mail addresses: iancu@dsm-mail.cea.fr (E. Iancu), mclerran@quark.phy.bnl.gov (L. McLerran). 
2 Membre du Centre National de la Recherche Scientifique (CNRS), France. 



Preprint submitted to Elsevier 



2 February 2008 



1 Introduction 



Recently, there has been significant progress in our understanding of the dynamics in pertur- 
bative QCD at high energy/high gluon density, leading to a consistent picture of the high-energy 
evolution as a classical stochastic process of a special type [1,2]: a non-local generalization of 
the reaction-diffusion process of statistical physics, with the 'non-locality' referring both to the 
transverse momenta and to the position in the two-dimensional 'impact-parameter space' — 
the plane transverse to the collision axis. However, the consequences of this new picture for the 
dynamics in impact-parameter space are still to be explored: although the impact-parameter 
dependence is in principle encoded in the underlying evolution equations — the Pomeron loop 
equations of Refs. [2-4] — , this dependence turns out to be too complicated to deal with in 
practice, and so far it has been neglected in the applications of these equations. Some important 
aspects of this dynamics — like the peripheral dynamics responsible for the Froissart growth of 
the total cross-section (see, e.g., [5,6]) — are clearly non-perturbative. Still, with increasing 
energy, there is an increasingly large region around the center of the hadron where the gluon 
density is high and perturbation should apply, including for a calculation of the correlations in 
b (the two-dimensional vector denoting the impact parameter). 

For long time, this high-density central region has been assumed to be quasi-homogeneous 
— a disk which appears uniformly black to any external probe whose resolving power Q 2 in 
the transverse space is smaller than some critical value fixed by the gluon density in the target. 
Such a uniform disk would be characterized by just two quantities: (i) the (local) saturation 
momentum Q 2 (b), which fixes the critical scale for 'blackness' (i.e., for the onset of unitarity 
corrections) and is roughly independent of b (it typically has a smooth profile decreasing from 
the center towards the periphery on a radial distance set by the 'soft' QCD scale Aq CD ), and 
(ii) the radius R of this central region where the density is high. Both quantities are expected to 
grow with the energy, but whereas the respective growth is rapid (power-like) for the saturation 
momentum and can be computed in perturbation theory — since determined by quasi-local 
gluon splitting processes within the high-density region — , that of the 'black disk' radius is 
much slower (logarithmic in s) and also non-perturbative, since it is related to the expansion of 
the black disk into the outer corona at relatively low density. 

Very recently, however, this traditional picture has been challenged [2, 7, 8] by the newly 
developed picture of the stochastic evolution, which suggests that, for sufficiently high energies 
at least, the central region at high density should be the site of wild fluctuations, leading to the 
coexistence of regions where the gluon density is anomalously high — i.e., the local saturation 
momentum Q 2 s (b) is considerably larger than the respective average value Q 2 — with regions 
which are anomalously dilute, such that Q 2 (b) <C Q 2 . The high-density regions appear as black 
spots to an external probe with a resolution Q 2 ~ Q 2 , whereas the low-density ones rather look 
'grey' on the same resolution scale, with the precise nuance of 'grey' depending upon the ratio 
between the local saturation scale Q 2 (b) and its average value Q 2 (see Fig.Q]). Remarkably, it 
turns out that, at least for a simple projectile like a dipole, the 'black spots' completely dominate 
the average scattering amplitude (T(Q 2 )) up to very high resolution scales Q 2 , well above 
Q 2 [2]. That is, even for values Q 2 which are so high that the scattering is weak on the average, 
i.e., such that (T(Q 2 )) <C 1, the average amplitude is still controlled by the dense fluctuations 
for which T(Q 2 ) ~ 1, i.e., by the rare events in which the projectile has hit a black spot. 

This new picture has rich and interesting consequences: it leads to a total breakdown of the 



2 




Fig. 1. A transverse view of a high-energy hadron, as 'seen' in deep inelastic scattering at very high 
energy. The dipole projectile probes a target area of the order of its own transverse area, at various 
impact parameters. 

'twist expansion' up to very high values for Q 2 and it predicts a new scaling law for (T(Q 2 )) 
at high energies [1, 2, 9] — known as diffusive scaling [7] — , which should replace the geo- 
metric scaling expected [10-12] from mean field approximations like the Balitsky-Kovchegov 
equation [13, 14] or the more general (functional) JIMWLK equation [15-17]. 

The strong fluctuations in the gluon distribution at high energy find their origin in gluon- 
number fluctuations in the early stages of the evolution, i.e., BFKL-like [18, 19] splitting pro- 
cesses^ which induce correlations between the gluons having a common ancestor, and thus 
between the 'spots' generated via the subsequent evolution of these gluons. So far, the effects 
of these correlations have been studied only at fixed impact parameter, via a coarse-graining of 
the Pomeron loop equations in impact-parameter space [1,2]. This approximation is perhaps 
sufficient for a study of the scattering of a small projectile which is quasi-localized in b (e.g., 
the case of the dipole scattering, as relevant for 7*p deep inelastic scattering [7], and also for pp 
or pA collisions under specific circumstances [8]), but on the other hand it prevents one from 
studying the correlations in b, which could be experimentally accessed via multiple-particle 
production in hadron-hadron collisions (say, at LHC). Moreover, as we shall later argue, this 
approximation has some shortcomings even for the description of the short-range correlations, 
for which it was a priori intended: by ignoring any information about the size of the fluctuations 
in impact-parameter space, it artificially suppresses the fluctuations with very small sizes, i.e., 
the tiny spots where the local saturation momentum is much larger than the average one. Or, 
as alluded to before, a proper description of such tiny spots is essential in order to study the 
scattering of small projectiles with resolution Q 2 ^> Q 2 . 

In this paper, we shall make a first attempt to go beyond the coarse-graining approximation 
introduced in Refs. [1,2]. Namely, we shall propose an effective field theory describing the dis- 
tribution of the saturation momentum in impact-parameter space at a fixed (high) energy and 
for a large number of colors N c . (The large-iV c approximation is needed in order to be able to 



3 The 'black spots' have been originally observed [20] in numerical simulations of Mueller's dipole 
picture [19] — the large- N c version of the BFKL evolution — , but at that time it was not clear whether 
they would survive after including the 'unitarity corrections' (i.e., the saturation effects due to the non- 
linear gluon dynamics at high density [13, 15, 17,21]). 
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neglect long-range color exchanges between saturated spots; see Sect. [2] for details.) Note that 
the saturation momentum Q s (b) is also the typical transverse momentum of the gluon configu- 
ration located around b. Hence, the effective theory that we shall propose encodes information 
about the gluon distributions in both momentum space and impact-parameter space. 

Let us start by emphasizing that the very existence of an effective field theory for Q s (b) is 
not at all obvious. For instance, it is not a priori clear whether one can effectively replace (even 
if only approximately) the information about the gluon distribution or the dipole scattering 
amplitudes by a theory for the distribution of the saturation momentum alone. Furthermore, 
even assuming that such a theory exists, it is not at all clear whether it can be made local. The 
Pomeron loop equations are non-local in the transverse momenta and coordinates, and this may 
translate into non-localities in b in the effective theory that we are looking for. Moreover, the 
correlations induced by the evolution are generally non-local in rapidity, since associated with 
gluon splittings in the intermediate steps of the evolution. Still, as we shall argue in Sect. |2l the 
short-range correlations, at least, are quasi-local in rapidity, since predominantly produced via 
gluon splittings in the late stages of the evolution, close to the final rapidity Y. By 'short-range' 
we mean distances of the order of, or smaller than, the average saturation length 1/Q (Y) 
at rapidity Y (with the notation^ Ql = (Q 2 ))- Hence, we expect a local field theory to be 
meaningful, at least, for the dynamics over such short distances. 

To proceed, we shall simply assume that a local field theory for Q s (b) exists and then we shall 
constrain its structure from general physical considerations. More precisely, we shall construct 
this theory as the natural generalization of the 6-independent distribution proposed in Refs. 
[1,2] which is consistent the uncertainty principle and the conformal symmetry of the high- 
energy evolution equations. Together, these constrains almost uniquely fix the structure of the 
effective theory, as we explain now: 

As already mentioned, the saturation momentum Q s (b) is also the typical transverse mo- 
mentum of the gluon configuration centered at b. By virtue of the uncertainty principle, the 
area occupied by that configuration cannot be smaller than \/Q 2 s {b). This area cannot be much 
larger either, since a gluon configuration which has reached saturation on some scale Q s (b) can 
hardly emit gluons with soft momenta k± <C Q s {b) [22-25]. That is, its subsequent evolution 
predominantly proceeds via the emission of harder gluons with k± Q s (b), which cannot in- 
crease the size of the configuration (except very slowly, via peripheral emissions). In turn, such 
hard gluons become spots which evolve towards saturation on their own sizes and at the same 
time act as sources for even smaller spots, which spread over the surface of the original config- 
uration and where the local saturation momenta are much harder. We see that the high-energy 
evolution is strongly biased towards smaller sizes, thus leading — in a three-dimensional pic- 
ture where the impact-parameter space occupies the (x, y)-plane and the saturation momentum 
Qs{b) is represented along the z-axis — to a 'landscape' picture, with spikes of various heights, 
randomly distributed and surrounded by valleys. 

Furthermore, the equations describing the evolution in QCD at high energy and in the lead- 
ing logarithmic approximation have the important property of conformal symmetry, that is, they 
are invariant under conformal (Mobius) transformations in the transverse plane. This has since 
long been appreciated in the context of the linear, BFKL, evolution [18], with important con- 



4 Note that so far we have introduced two different notations for the 'average saturation momentum', 
namely Q 2 S and Qq = (Q 2 .); they correspond to different definitions, as we shalf expfain in Sect. 2. 
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sequences [26-28] (in particular, for the calculability of the theory; see Ref. [29] for a review 
and more references), but this remains true for the non-linear BK [13, 14], JIMWLK [15-17], 
or Pomeron loop [2-4], equations, since all these equations involve the BFKL splitting kernel 
(for gluons or dipoles) together with gluon-number changing vertices, which indeed respect 
conformal symmetry. One may think that this symmetry is inconsistent with the emergence of 
a special scale — the (average) saturation momentum — in the solutions to these equations, 
but this is actually not true: the saturation scale is not spontaneously generated by the evolu- 
tion, rather it comes out as the evolution of the scale introduced by the initial conditions at 
low energy [10, 11,21,23]. On the other hand, the gluon correlations are generated by the 
conformally-invariant evolutioiO], so we expect this symmetry to constrain the effective theory 
describing these correlations. 

The effective theory that we shall arrive at by exploiting such considerations is a two- 
dimensional, interacting, scalar field theory, in which the field <\> is proportional to the logarithm 
of the local saturation momentum, and which involves two 'free' parameters: the expectation 
value Ql = (Q 2 ) of the saturation momentum and a dimensionless coupling constant a, which 
characterizes the degree of disorder introduced by gluon-number fluctuations in the course of 
the evolution. Both parameters are expected to rise with the energy, in a way which is in prin- 
ciple determined by the underlying QCD evolution equations, but which needs not be specified 
for our present purposes. It suffices to say that the weak coupling regime a <C 1 in the effective 
theory corresponds to low, or intermediate, energies, where the dynamics is quasi-deterministic, 
whereas the strong coupling regime o ^> 1 corresponds to the more interesting situation at high 
energy, where the stochastic aspects are fully developed and essential. 

The effective theory involves three basic ingredients: (i) the standard kinetic term, which 
couples fluctuations at neighboring points, (ii) a potential term, which ensures that the typical 
gradients at b are of order Q 2 (b), in fulfillment of the uncertainty principle, and (iii) a source 
term, which enforces the value Ql = (Qg(b)) of the average saturation momentum. Since 
proportional to Q 2 S oc e "*, the potential is exponential in 0, which is precisely as it should 
for consistency with conformal symmetry: together, the kinetic plus the potential terms are 
recognized as the Liouville action, which is conformally-invariant, and thus integrable. The 
quantum Liouville field theory (LFT) has been extensively studied over the last decades — 
especially, in connection with studies of quantum gravity in two dimensions — and many exact 
results are known by now about its properties [30-35]. However, precisely by virtue of its 
symmetry, LFT involves no mass scale (it gives rise to power-law correlations on all distance 
scales), and thus cannot accommodate a non-zero value for the average saturation momentum. 
This is taken care off by the source term, linear in cf>, at the expense of explicitly breaking the 
conformal symmetry. This breaking has dramatic consequences on the properties of the theory 
and it also complicates its analysis very much (since it is not possible to directly exploit the 
large amount of information known about LFT). 

Yet, some general properties of the effective theory can be inferred by semi-classical and 
perturbative techniques, or simply by inspection of the action. The presence of the source term 
endows this theory with a stable ground state (in contrast to LFT), which appears as a saddle 
point of the action and allows for a perturbative treatment of the weak coupling (a C 1) regime. 
The perturbative calculations, that we shall push up to two-loop order, reveal a very interest- 



5 We assume that there were no correlations in the initial conditions at low energy, for simplicity. 
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ing property, which is furthermore confirmed, at non-perturbative level, by inspection of the 
corresponding Dyson equations: the correlations of the saturation momentum (an exponential 
operator in this theory) are ultraviolet finite, and not just renormalizable, as generally expected 
for a two-dimensional field theory. This means that the operator Q 2 S has no anomalous dimen- 
sion, which in turn implies that its short-range correlations, over distances R < 1/Qo, have a 
power-like decay, with exactly the same powers as in pure Liouville theory — since these pow- 
ers are fixed by the natural dimension of the operator. On the other hand, the theory predicts 
that these correlations decay exponentially on larger distances R ^> 1/Qo, with the (average) 
saturation momentum Q playing the role of a screening mass. The emergence of an exponen- 
tial fall-off in the context of perturbative QCD may look surprising, but it is presumably an 
artifact of our insistence on a field theory which is local in rapidity: the physical correlations on 
large distances R^> 1/Qo are typically generated via gluon splittings in the early stages of the 
evolution (the earlier, the larger R is; see Sect. [2])), which are not encoded in our present theory. 

The paper is organized as follows: Sect. [2]presents a brief and critical discussion of the b- 
independent distribution proposed in Refs. [1,2], with the purpose of clarifying its limitations 
and, more generally, the limitations of any effective theory which is local in rapidity. In Sect. 
[3] we construct the low-energy/weak coupling (a < 1) version of the effective field theory, 
as the straightforward extension of the 6-independent distribution alluded to above. Sect. |4] 
is our main section: after a quick introduction to the Liouville field theory and its conformal 
symmetry, we explain the relevance of this theory for the QCD problem at hand, then present 
the complete action for our effective theory (including the symmetry-breaking source term) and 
discuss some of its properties. Sect. [5] is slightly more technical, since devoted to perturbative 
calculations to two-loop order, with the purpose of demonstrating the ultraviolet-finiteness of 
the correlation functions when computed in the effective theory. Finally, Sect. 6 summarizes 
our results and conclusions. In the Appendix, we show how to couple the effective theory to 
the CGC formalism (within the simple context of McLerran-Venugopalan model [36]), which 
is useful in view of computing observables. 

2 The coarse-graining approximation and its limitations 

Our starting point is the Gaussian probability distribution for the saturation momentum intro- 
duced in Refs. [1,2] (see also Refs. [37-39]), which can be viewed as a coarse-grained version 
of the effective theory that we intend to construct, in which the impact-parameter dependence 
has been averaged out. The applicability of such a coarse-graining will be shortly discussed. 

The random variable in this distribution is the logarithm p s = \n(Q 2 / A 2 ) of the saturation 
momentum (with A some arbitrary scale of reference), which is the scale which separates, in an 
event-by-event description, between a high-density, or 'Color Glass Condensate' [40], phase at 
low transverse momenta (k± < Q s or p = ln(/c^/A 2 ) < p s ), where the gluon occupation factor 
saturates^ at a large value of 0(l/a s ), and a low density phase at high momenta (p > p s ), 
where the gluon occupancy is still low but it rises rapidly with the energy (as a power of s), via 
BFKL gluon splitting. It is customary to work with the 'rapidity' variable Y = In s, which plays 
the role of an 'evolution time' for the high-energy evolution. When increasing Y, modes with 

6 Strictly speaking, the gluon occupancy keeps growing with the energy in the CGC phase, but only very 
slowly: logarithmically in s, that is, linearly in the rapidity Y = Ins [23,24]. 
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Fig. 2. A picture of the diffusive saturation boundary in the presence of fluctuations; shown are the 
average saturation line (p s ) = Xa s Y, the diffusive radius a oc VY for the Gaussian distribution (12.11) . 
and the diffusive scaling window with a width a 2 oc Y. 



higher and higher values of kj_ enter at saturation, so the borderline p s propagates towards higher 
transverse momenta. However, due to gluon-number fluctuations in the splitting process — 
which can be associated with the fact that the particle number is discrete [1] — , this progression 
of p s is not uniform but stochastic, like a one-dimensional Brownian motion. 

Based on the analogy with the reaction-diffusion problem [37-39,41], it has been argued in 
Refs. [1,2] that p s should be a Gaussian random variable with probability distribution 
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(2.1) 



where both the central value (p s ) and the dispersion a 2 = (p 2 ) — (p s ) 2 are expectecCO to rise 
linearly with Y (for Y large enough): (p s ) = \a s Y and a 2 ~ Da s Y, with a s = a s N c /n. The 
coefficients A and D are known [37,38] only in the formal limit a s — > 0, and the corresponding 
expressions cannot be extrapolated to realistic values of a s since they depend upon ln(l/a^). 
For what follows, their actual values are unimportant, and so are the precise F-dependencies 
of (p s ) and a; all that matters is that these quantities rise quite fast with Y. So long as a < 1 
— the case for relatively low energies — the theory remains quasi-deterministic. But the most 
interesting regime for us here is the high-energy regime at a > 1, where the fluctuations are 
fully developed and have a strong influence on the theory. 

Since the saturation momentum is fluctuating, a physical discussion will naturally involve a 



7 Recall that, throughout this analysis, we consider the leading-order approximation where the coupling 
is not running. 
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range of values for Q s , rather than just a single value (see Fig. [2]). The range which is a priori 
privileged by the Gaussian probability (12.11 ) is the 'diffusive disk' at \p s — (p s )\ < a, where 
the probability is reasonably large: P(p s ) ^ \ jo. But the relevant range also depends upon the 
physical quantity of interest, which can favor some values of Q s over the others. The simplest 
example in that sense refers to the very definition of the 'average saturation momentum', which 
is not unique (because the relation between Q s and p s is non-linear: Q 2 = A 2 e Pa ) and thus 
needs to be properly specified — different definitions can be relevant for different problems. 
The form of the Gaussian distribution (12.11 ) makes it natural to define 

Q2 = A 2 e (p s ) _ (2 2) 

But the expectation value of the 'operator' e Pa can be also computed, with the following result: 

(Ql) = A 2 (e pa ) = A 2 e {ps) e^ 2 = Q 2 S ^ . (2.3) 

This is always harder than the scale in Eq. (12.21) . and for sufficiently high energies, where a ^> 1, 
it is even much harder. This considerable difference between the two definitions can be easily 
traced back: in computing (Q 2 ), the exponential operator e pa biases the integration towards 
very large values of p s , of order p s ~ (p s ) + a 2 , as opposed to the typical value p s ~ (p s ) 
contributing to Q 2 . Note that the large deviation p s — (p s ) ~ a 2 corresponds to fluctuations 
which are very rare: P(p s ) ~ exp(— a 2 /2) <C 1, and which occupy a very small area ~ 1/ (Q 2 ) 
in impact-parameter space (by the uncertainty principle). Yet, such rare and tiny fluctuations 
are physically important, as they define the upper bound of the diffusive scaling window [2,7,8]: 

2 2 

- y « P - (Ps) « y , (2.4) 

(see also Fig. [2]) where p = \n(Q 2 /A 2 ), with Q 2 the transverse resolution of a small projectile, 
like a 'color dipole' (a gg-pair fluctuation of the virtual photon in deep inelastic electron-hadron 



scattering), which scatters off the hadron. The interval in Eq. (12.41) represents the kinemati- 
cal window at high energy (cr 3> 1) within which the (average) dipole scattering amplitude 
(T(p,Y)) 'scales' as a function of the dimensionless variable Z = (p — (p s ))/a rather than 
separately depending upon p and Y : (T(p, Y)} w (T(Z)) (see Refs. [2,7, 8] for details). The 
diffusive scaling window in Eq. (12.41) can be alternatively described as 



Qie-v « Q z « Qie>° ee (Qi), (2.5) 

which explicitly shows that both definitions for the 'average saturation momentum' introduced 
in Eqs. (|2.2I) and, respectively, (12.31) play a role in characterizing the physical effects of the 
fluctuations in Q s . 

At this point one should remember that the previous discussion was based on the local ap- 
proximation (12.11 ) for the distribution of Q s , so it is important to understand what are the as- 
sumptions beyond this approximation. In Refs. [1,2] (see, especially, the discussion in Sect. 6 
of Ref. [2]), the distribution (|2.1I) has been obtained via a coarse-graining in impact-parameter 
space: the b-dependence of the dipole scattering amplitudes, as encoded in the Pomeron loop 
equations [3,4], has been averaged out over a region in the impact-parameter space which in 
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Refs. [1,2] has been loosely characterized as the 'dipole size', but which should be more prop- 
erly interpreted as an intrinsic scale in the hadron, so like 1/Q S . In this averaging, one has 
potentially neglected two types of correlations: (i) short-range correlations between the fluc- 
tuations ('spots') with small sizes R <C 1/Q S (i.e., with large saturation momenta Q s 3> Q s ) 
which lie inside the coarse-graining cell — such smaller spots were treated as being uniformly 
distributed over the cell area l/Q 2 , and ( i i ) the long-range (R 3> 1/ Q s ) correlations between 
different cells. Besides, one has assumed homogeneity in b — the average saturation momen- 
tum Q1 and the dispersion a 2 were taken to be independent of b — , which is not really an 
approximation, but merely a choice for the initial conditions at low energy (the homogeneity 
being preserved by the evolution according to the Pomeron loop equations). Under these as- 
sumptions, the 6-dependence has disappeared from the evolution equations, which were then 
shown [2] to be equivalent to a stochastic equation of the sFKPP type — the Langevin equation 
for the reaction-diffusion process [41]. The Gaussian probability distribution in Eq. (12.11) then 
follows from known properties of the sFKPP equation, as recently clarified in the respective 
literature [37-39]. 

Thus, the fact that here is no explicit b-dependence in Eq. (12.11) must be understood as fol- 
lows: this formula is meant to apply to points b within a given coarse-graining cell, with an area 
~ 1/Os '■> at ai l sucn points, the event-by-even local saturation scale is assumed to be the same, 
and equal to p s (in logarithmic units); that is, all the points within a cell fluctuate coherently 
with each other, so like a rigid body. On the other, nothing is said about the correlations be- 
tween different cells: Eq. (12.11) apply only to local fluctuations within a given cell, as averaged 
over the size of that cell. 

How does this approximate picture compare with the physical reality ? There are a pri- 
ori two mechanisms for building spot-spot correlations: (a) Long-range interactions between 
saturated spots: A gluon configuration with saturation momentum Q s has a low but non-zero 
probability to emit low-momentum gluons with k± <C Q s , thus creating long-range color fields 
over distances R ^> 1/Q S . These are dipolar fields, and not Coulomb ones, because the total 
color charge gets screened at saturation [24,25,42]. These fields mediate dipole-dipole interac- 
tions between saturated spots which are far away from each other. However, these interactions 
are suppressed by a factor l/N 2 , as they imply color exchanges, and thus can be safely ne- 
glected in the large-iV c approximation underlying the present discussion, (b) Gluon-number 
fluctuations in the high-energy evolution: As mentioned in the Introduction, different spots can 
be correlated with each other because they have a common ancestor at some earlier rapidity. 
Such correlations survive at large N c , since associated with fluctuations in a colorless quantity: 
the number of gluons in the light cone gauge. Eq. (|2.1I) is meant to capture these fluctuations 
in the coarse-graining approximation, and our purpose here is to relax this approximation by 
restoring the b-dependence of the induced correlations. 

However, as also mentioned in the Introduction, the relevant correlations are generally non- 
local in rapidity, since induced through branching processes which can occur at all the interme- 
diate rapidities. As we shall explain no the short-range correlations generated in this way 
can nevertheless be encoded in a local (in Y) field theory, of the type we would like to construct. 
To that aim, we shall consider two spots which are separated by a distance R at the final rapidity 
Y . These spots are correlated with each other provided they had a common ancestor at some 



We would like to thank Al Mueller for illuminating discussions on this particular issue. 
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earlier rapidity < y < Y. It turns out that y is correlated with R. Indeed, the typical spots at 
rapidity y have a saturation momentum Q s (y) and a size ~ 1/Q s (y). If that size is much smaller 
than R, then such a typical spot has only little probability to emit a gluon at a distance R away 
from its center and thus initiate the evolutions leading to the two final spots that we measure 
at Y. (This probability decreases as a large power of (RQsiy))^ 1 because of the screening of 
the color charge at saturation.) If, on the other hand, 1/Q s (y) ^> R, then there is a geometrical 
penalty factor ~ (RQ s (y)) 2 for the gluon to be emitted inside the (relatively small) domain 
with area R 2 . Hence, the most important intermediate rapidity y for creating correlations over 
a distance R is the one for which Q s (y) ~ 1/R. Therefore, the larger is R, the more we have 
to go backwards in the evolution, and the more non-local are the respective correlations in Y. 
Vice-versa, the short-range correlations with R ^ 1/Q S (Y) are typically produced in the late 
stages of the evolution, at y ~ Y, and hence are quasi-local in rapidity. It should be therefore 
possible to encode such short-range correlations in a local, effective, field theory. 

As we shall later see, these considerations are indeed consistent with the effective field theory 
that we shall arrive at, except for the fact that the actual scale for correlations which will emerge 
from that theory is not the scale Q 2 that we have focused on in the above discussion, but rather 
the harder scale (Q 2 S ), cf. Eq. (1231) . 



3 The effective field theory: weak coupling 

With this section, we start our program aiming at extending the coarse-grained distribution 
in Eq. (12.11) to an effective field theory which describes the distribution of the saturation mo- 
mentum in impact-parameter space. For more clarity, we shall first develop our arguments for 
the low energy/weak coupling regime d < 1, where we shall argue that the corresponding ex- 
tension is a free field theory for a massive scalar field in two (Euclidean) dimensions. Then, in 
Sect.Hl we shall present the generalization of this theory to arbitrary values of the coupling o. 

For more clarity, we shall assume 'mean-field-like' initial conditions at low energy (Y = 0): 
the initial gluon density is large and homogeneous, and the saturation momentum Q S (Y = 0) 
takes the same value A (with A 2 ^> Aq CD ) at all the points. By boosting this system to high 
energy, the gluon distribution evolves via (generally, non-linear) gluon splitting, and inhomo- 
geneities ('spots') appear in the event-by-event description, due to gluon-number fluctuations. 

In the early stages of this evolution — namely, so long as a <C 1 — , the fluctuations have no 
time to significantly develop, so all the spots have more or less the same siz e0,~ 1/Q s ,andthe 
same value for the saturation momentum, ~ Q s ; furthermore, they are only weakly correlated 
with each other. It is then straightforward to extend the coarse-grained approximation (12.11) to a 
distribution which covers the whole transverse profile of the hadron: this is simply the product 
of independent Gaussian distributions like that in Eq. (12.11 ), one for each spot : 



py[ps\ = n-7=- ex p 

„• v In a 



(Ps{j) - (Ps)) 2 



for o-<l, (3.1) 



2a 2 

where the discrete variable j labels the spots, and each spot has roughly an area l/Q 2 . Intro- 



9 Note that, when a <C 1, the various 'average saturation scales' characterizing the statistical ensemble 
(cf. Sect. 12) are close to each other: Q 2 m (Q 2 ). 
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ducing the fluctuation field r/(j) = p s (J) — (p s ), Eq. (13.11) implies: (rj(j)rj(l)) = a 2 5ij. So far, 
the quantities p s (j) or rj(j) refer globally to a spot. In order to promote them to field variables 
p s (x) or T](x), with x denoting the impact parameter, but keep the same correlations as above, 
one needs to introduce a kinetic term to smear out correlations over the spot size ~ l/Q s - The 
simplest kinetic term, and also the only one to be renormalizable in two dimensions, is the 
standard, quadratic, kinetic term, that we shall adopt in what follows. We are thus led to the 
following field-theoretical extension of Eq. (13.11 ) 

P Y [ V } = exp {-So{v]} , with So{v] = ^2 J d " x {(^vf + QW} , (3-2) 

which provides the weight function for the functional integral over 77(33) which defines expecta- 
tion values in the framework of this effective field theory. The measure in the functional integral 
is assumed to carry the proper normalization: / D[q] exp {— So [77]} = 1. 

Eq. (13.21 ) implies (r)(x)r](y)) = a 2 D(x — y), with D the propagator of a free, massive, 
scalar field in a two-dimensional Euclidean space: 

where r = \x — y\ and K is the respective Bessel function, with the following limiting be- 
haviours at short and, respective, large distances: 



K Q (z) 



In - , for z <C 1 



z 



(3.4) 



I V 2z 



^-e \ for z > 1. 



We see that the correlations generated by the field theory in Eq. (13.21) have indeed the sought for 
structure: they are quasi-uniform on relatively short distances, r ^ l/Q s , (i.e., among points 
located within a same spot), but they decay very fast (exponentially) over distances much larger 
than the typical size of a spot. 

The emergence of exponentially decaying correlations at large distances should be taken 
with a grain of salt (cf. the discussion in the Introduction): It corresponds to the fact that the 
correlations over larger distances r ^> 1/Q S are predominantly produced via splittings occuring 
in the early stages of the evolutions, i.e., at rapidities y < Y, whose effects are not included in 
the present formalism, local in Y. As for the late splittings (those taking place at y ~ Y), they 
produce correlations which fall off according to power laws at large separations r 3> 1/Q S , 
because the emissions of soft gluons with k± <C Q S (Y) is strongly suppressed by saturation. 
The fact that, in our formalism, this fall-off appears to be exponential, rather than power-like, 
is because we have not provided a faithful description of the gluon spectrum, including its 
softening at momenta below Q s , but rather we have replaced this whole spectrum by a unique 
scale — the saturation momentum — , which represents the average transverse momentum of the 
saturated gluons. To summarize this argument, the theory with action (13.21) provides the correct, 
logarithmic, behaviour at short distances, and it mimics the rapid fall-off at larger distances by 
an exponential tail, rather than the correct, power-law, one. This will be a generic feature of the 
effective theory, including at strong coupling. 
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We now return to the functional distribution in Eq. (13.21) and consider its predictions for the 
correlation functions of the operator saturation momentum, Q 2 (x) = A 2 e pa = Q 2 s e v . One finds 

(Q 2 (x)) = Q 2 e*^°) , 
{Ql{x)Ql{y)) = Qt e^ +D ^ , (3.5) 

etc. Strictly speaking, these expressions should be trusted only up to terms of 0(<j 2 ), since 
we are working under the assumption that cr C 1; notwithstanding, we here display the full 
exponentials, for comparison with subsequent results at strong coupling. 

The above formulae involve the equal-point limit -D(O) of the propagator, which is diver- 
gent: -D(O) = ( l/27r) ]n(l/aQ s ), with a a short-distance cutoff ('lattice spacing') introduced 
to regularize the divergence. We see that, even within this free field theory, the exponential (or 
'vertex') operator V(x) = exp{i](x)} develops ultraviolet divergences, which arise from con- 
tractions internal to this operator. As is well known [43], and also manifest on Eq. (13.51 ), such 
divergences can be eliminated via multiplicative renormalization of the vertex operator (which 
here is tantamount to normal-ordering the polynomial operators which appear in its expansion). 
Namely, we shall define the renormalized vertex operator as 

V R (x) = V(x)exp | ^ In (ae 2 *Q s ) j = exp [1 - D(0)]| , (3.6) 

where the particular subtraction point 1/a = e 2n Q s has been chosen for reasons to shortly 
become clear. One then obtains the following, finite, vertex correlation functions: 

(V R (x)) = e^ 2 , ... , (V R ( Xl )---V R (x n )) = expjV + a 2 ^/^-^)}, (3.7) 



which in turn imply, for the renormalized operator Q 2 (x) = Q 2 Vr(x), 
(Q 2 s (x)) = Q 2 S e^ 2 , 

(Q 2 s (x)Q 2 s (y)) = Qt e^ 2D ^ = (Q 2 S ) 2 exp{a 2 D(x - y)} . (3.8) 

One can now appreciate the particular choice for the subtraction point in Eq. (13.61) : this is such 
that the expectation value (Q 2 ) computed in the present field theory precisely matches the cor- 
responding prediction of the coarse-grained approximation, cf. Eq. (12.31) . 

Interestingly, Eq. (|3.8I) reveals the emergence of power-like correlations over (relatively) 
short distances, that is, in between the points lying inside a same spot: 

(O 2 ) 2 

(Q 2 s (x)Q 2 s (y)) « Ws/ _ ict2/2 ^ for \x - y\ « 1/Q S . (3.9) 

\x - y\Q s 

This exhibits a singularity in the equal-point limit, which is however physical, in the sense that 
one cannot localize fluctuations in a quantum field theory down to a point without generating 
singularities. On the other hand, at large distances, the connected piece of the 2-point function 
in Eq. (13.81) dies away exponentially: 
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(Q 2 s (x)Q 2 s (y)) - (Q 2 ) 2 oc a 2 (Q 2 ) 2 exp(-|x - y\Q s ) for \x-y\^>l/Q s . (3.10) 



Once again, this exponential decay should be taken with a grain of salt: the effective theory is 
not supposed to apply to such large separations. 

4 The effective field theory: general case 

Although this has not been explicitly spelled out in the previous discussion, the free field 
theory of Eq. (13.21) is indeed consistent with the uncertainty principle in the regime where this 
theory is meant to apply (i.e., for a <C 1): indeed, the kinetic term there smears out inhomo- 
geneities over the very short distances r <C 1/Q s ,'m agreement with the fact that the size of a 
spot cannot be smaller than the inverse of the typical momentum of the gluons composing that 
spot, namely kj_ ~ Q s - However, for sufficiently high energy, such that a > 1, the fluctuations 
become important and then the actual saturation momentum Q 2 (x) at a given point and in a 
given event can be very different from any of its expectation values {Q 2 or (Q 2 )) previously 
introduced. In such a case, the uncertainty principle requires that the minimal size of a spot 
located at x be fixed by the actual value Q 2 (x) of the saturation momentum, and not by its 
expectation value. This means that, in the general action for r/, the kinetic term should compete 
with the actual (event-by-event) saturation momentum, and not with its 'expectation value' 
(whatever the meaning of the latter is). 

Accordingly, the effective 'mass term' in the action should be the field-dependent scale 
Q 2 (x) = Q 2 exp{r](x)}. This argument suggests the following generalization of Eq. ( 13.21 ): 



which involves an exponential potential. Clearly, this is now an interacting field theory in which 
the parameter a — which, we recall, is a measure of the dispersion introduced by fluctuations 
— plays the role of a coupling constant, as it can be recognized after rescaling 77 = a<p '■ 



Strong interactions for the field correspond to strong correlations between spots with different 
sizes and at different locations, as physically expected at sufficiently high energy. However, it 
should be clear that the precise form of the potential for (ft, beyond the exponential factor, cannot 
be uniquely fixed by the uncertainty principle alone: a factor of (j) counts like In Q s and hence 
it cannot modify the power-counting argument that V 2 ~ Q 2 . Thus, the uncertainty principle 
alone would allow for any potential where e ^ is multiplied by an arbitrary polynomial in </>. 

At this point we shall invoke conformal symmetry to further specify the form of the potential. 
As mentioned in the Introduction, the high-energy evolution equations in QCD in the leading 
logarithmic approximation are invariant under scale, and, more generally, (special) conformal 
transformations. It is likely that a similar symmetry should also hold for the correlations gen- 
erated by this evolution, at least within limited ranges. This is a very strong constraint on the 
effective theory, which almost uniquely fixes its structure, as we shall explain in what follows. 




(4.1) 




(4.2) 
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4.1 Liouville field theory in a nutshell 



Let us temporarily assume that the effective theory should have exact conformal symmetry 
(this assumption cannot right, as we shall later argue, but it allows us to provide a first iteration 
for the effective action). Then, the potential for is necessarily a pure exponential, and the 
corresponding effective action is the same as the classical Liouville action 10 1 [31,35] : 

S L [4>] = J d 2 x |^(V» 2 + 3f e*J . (4.3) 

A factor l/cr 2 has been introduced in the coefficient of the potential to ensure that, in the weak- 
coupling limit a — ► with fixed r\ = a<p, both terms in the action are of 0(1/ a 2 ). Also, for 
later convenience, we have replaced the mass scale Q 2 S in front of the potential by the generic 
scale Qq. Correspondingly, we re-interpret the field <fi as 

0(a) = I In (4.4) 

so that the quantity Q^e* 7 ^ 33 ) = Q 2 s (x) preserves its original meaning as the event-by-event 
saturation momentum (a composite, 'vertex', operator in the present field theory). 

The exponential potential has the remarkable property to be invariant under conformal trans- 
formations provided one allows the field <\> to change by a shift under such transformations. 
Consider indeed the scale transformation x — > Xx with A > 0. Then the action (14.31) is invariant 
under the following transformations: 

x — > x = Xx 

<j>{x) -> 4>'(x') = <j>(x) - - In A. (4.5) 

a 

More generally, Sl is invariant under local conformal transformations, that is, the general holo- 
morphic transformations of the complex plane (which include the special, or global, conformal 
transformations). To formulate this symmetry, it is convenient to introduce complex notations: 
x = (x±, X2) — > z = xx + 1X2. Then, the action (14.31) is invariant under 



z -* w = f(z) 

(z) -> <j>'(w) = <f>(z) - - In 

a 



df 



dz 



2 



(4.6) 



since, e.g., e^dzdz = e^'dwdw. 

Via an appropriate quantization procedure, which turns out to be quite non-trivial [31, 35], 
the above symmetry properties can be carried over to the quantum version of the Liouville field 
theory (QLFT), with important consequences: First, being invariant under a continuum class 
of symmetry transformations, the quantum Liouville theory is integrable. Second, the point- 
dependence of the correlation functions of the vertex operator V(x) = exp{a(p(x)} is, to a 



10 The classical Liouville theory was extensively studied at the end of the nineteenth century in connec- 
tion with the uniformization problem for Riemann surfaces (see, e.g., the discussion in [32]). 
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large extent, fixed by the Ward identities associated with conformal transformations (see, e.g., 
the textbook discussion in [44]). This yields, e.g., 



(V(x)V(y))= i _ C *W , 



x - y\ A Qt 

C 3 (a) 

x ~ y\ 2 \y ~ z\ 2 \z — x 2 \Qq 



(V(xMy)V(z)) = - — > M , (4.7) 



where the coefficients C^er), C^{a), etc., are known exactly [33,34]. These formulae follow 
from the conformal symmetry of the (properly defined) path integral together with the fact that 
V{x) is a primary field with scaling dimension A = 1 (cf. Eq. (|4.5I) ) : 

x _> x > = Xx => V(x) -> V'(x') = V(x) = ^ V(x) . (4.8) 

This particular transformation law for the vertex operator is natural from the point of view of 
QCD: under a scale transformation, the saturation momentum Q 2 s (x) = QlV(x) transforms as 
an operator with mass dimension two, which is its physical dimension indeed. 

A rather subtle aspect of LFT, which complicates its quantum implementation, refers to the 
interplay between ultraviolet renormalization and conformal symmetry: although superrenor- 
malizable (since the exponential potential can be expanded out in a series in powers of <p and 
each term in this series is superrenormalizable in d = 2), the quantum theory requires renor- 
malization to remove 'tadpole' divergences — i.e., divergences associated with the equal-point 
limit of the propagator. Such divergences arise from contracting fields at the same point and thus 
can be eliminated by normal-ordering the operators. But this procedure introduces 'anomalous 
dimensions' for the vertex operators (their scaling dimensions acquire quantum corrections), 
which could spoil conformal symmetry. To maintain the symmetry, the theory is modified in 
such a way to ensure that the vertex operator V = which enters the action preserves at quan- 
tum level the classical dimension A = 1, cf. Eq. (14.81) . Note however that for a generic vertex 
operator V a = e°"^ with a ^ a, the ensuing quantum dimension is different from the respective 
prediction A^ = a/a of the classical Liouville theory (cf. Eq. (14.51) ) (see, e.g., [31,34,35]). 
Here, we do not need to discuss these complications in more detail because, as we shall see, 
they do not show up in the modified version of the Liouville theory that we shall propose as an 
effective theory for the QCD problem at hand. In fact, the only reason for us to mention these 
subtleties here, it is to emphasize, by contrast, the situation in our final theory (cf. Sect. 14.21 ), 
where all such UV complications disappear. 

But before turning to that presentation, let us rapidly explain why the standard QLFT cannot 
be, by itself, the complete effective theory that we need in QCD. The problem comes precisely 
from that feature of this theory which is also its main virtue : the exact conformal symme- 
try. Because of this symmetry, QLFT involves no real mass scale — the scale Ql apparent in 
Eq. (14.31 ) has no intrinsic meaning since its magnitude can be changed at will (and thus made 
arbitrarily large or arbitrarily small) by shifting the field </> under the path integral — , and hence 
it cannot describe a non-trivial gluon distribution, as characterized by a non-zero expectation 
value for the saturation momentum. In fact, the Ward identities for conformal symmetry [44] 
imply (V(x)) = 0, showing that, from the perspective of QCD, the LFT would be a theory for 
fluctuations, but ... without matter: (Q 2 s (x)) = at any x in pure QLFT ! 
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The same basic problem can be seen under different angles, revealing as many 'paradoxes' 
from the point of view of QCD: The exponential potential in Eq. (14.31) has no local minimum 
and becomes flat when <\> is negative and large; hence, the classical field is rolling down to 
— oo and, correspondingly, the quantum theory has no stable ground state. Because of that, it 
makes no sense to compute the correlation functions of the field <p — only the correlations 
of the vertex operator V a = (with generic a), or the mixed correlations involving both 
and V a , are a priori well defined. Still because of the lack of a stable ground state, there is no 
fundamental difference between 'weak' and 'strong' coupling in QLFT: the vertex correlation 
functions exhibit the same power-law behaviour, cf. Eq. (14.71) . for any value of a, small or large, 
since this behaviour is fixed by conformal symmetry alone. Besides, the coefficients in these 
correlations, so like C^er) and 03(a) in Eq. (|4.7I) . show interesting 'self-duality' properties 
under the exchange a <-> 1/a [34]. There is furthermore no distinction between 'short' and 
'long' distances, precisely because there is no intrinsic mass scale in the theory: QLFT generates 
power-law correlations for the vertex operators on all scales. 

Clearly, all these properties would be unacceptable in the framework of our original QCD 
problem, where, on the contrary, we expect pronounced differences between the weak and the 
strong coupling regimes, or between short and large distances. To understand the way out of 
such paradoxes, let us remind that, in the context of QCD, the conformal symmetry character- 
istic of the evolution equations is explicitly broken by the initial condition at Y = 0, which 
introduces a physical scale in the problem. It is the subsequent evolution of this scale with in- 
creasing energy which fixes the average saturation momentum at some later Y. But the effective 
theory that we are looking for is not a theory for the evolution, but rather for the results in this 
evolution (in terms of correlations of Q s ) at the final rapidity Y. In this theory, the average 
saturation momentum at Y is a parameter that must be introduced by hand, with the effect that 
conformal symmetry is explicitly broken. A particularly simple implementation of this idea will 
be described in the next section. 

4.2 Liouville field theory with a homogeneous source term 

From now on, we shall adopt the point of view that the Liouville action describes the fluc- 
tuations in the high-energy evolution (since it has the correct symmetry in that sense), but that 
in order to also describe the average gluon distribution, this action must be supplemented with 
a 'source term' which breaks down the scale symmetry and thus provides a non-zero expec- 
tation value for the saturation momentum. We shall consider here the simplest source term, 
namely, an operator linear in whose strength (the source density) is adjusted in such a way 
to produce a prescribed value for the average saturation momentum. For simplicity, we shall 
first consider a homogeneous gluon distribution, as produced by the high-energy evolution of 
an initial distribution which was itself homogeneous (say, a 'large nucleus' in the context of 
the McLerran-Venugopalan model [36]). This is not a very strong restriction, since the typical 
length scales that we shall be interested in at high energy are anyway much shorter than the 
scales characterizing the inhomogeneity in a physical hadron at low energy. 

We thus introduce the following extension of the Liouville action: S = Sl + j f x 4>(x), 
where Sl is given by Eq. (14.31) and j is the source density, assumed to be homogeneous. There 
are several ways to fix j, all of them leading to the same result. For instance, by requiring S to 
reduce to the Gaussian action (13.21 ) in the limit of small fluctuations/weak coupling a<p <C 1, 
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one immediately finds that the source term should remove the linear term from the expansion 
of the Liouville exponential. One thus finds: 

S[0] = / d 2 * |i(V» 2 + 5f (e°+ - o<t> - l) } , (4.9) 

where we have also subtracted the constant, unit, term from the exponential, for convenience. 
By construction, the above action has the saddle point = 0, which controls the dynamics in 
the weak coupling (i.e., low energy) regime a -C 1. More generally, the above potential has a 
unique minimum at = 0, which implies that the quantum theory defined by S has a stable 
ground state (at least, in perturbation theory). In particular, within this theory it makes sense to 
perturbatively compute the correlation functions of (unlike in LFT). 

The comparison between Eq. (14.91) and the quadratic action in Eq. (13.21) seems to suggest 
Ql = Q 2 S . However, this identification is wrong in general — it holds only as an approximate 
equality in the weak-coupling regime where the action (14.91) is supposed to reduce to Eq. (|3.2I) . 
To make the proper identification in the general case, consider the first Dyson equation gener- 
ated by Eq. (14T91) . that is 

{m)=° ~ «•"» 

By homogeneity, the mean field (4>(x)) is independent of x, hence V 2 (0) = 0, and then 
Eq. (14.101 ) implies: 

(e*) = l, or (Q 2 S ) = Ql (e*) = Q 2 . (4.11) 

Thus, the scale Ql in the action (14.91) must be understood as the average saturation momentum 
in the sense of Eq. (12.31) (and not of Eq. (12.21) !). In fact, one could have alternatively determined 
the value of the source density j by directly requiring (Q 2 ) = Ql (rather than going through a 
weak coupling argument); via the first Dyson equation, this condition would have again implied 
j = —Ql/a, as above. 

Note the importance of the source term in the action for the previous argument: without 
that term, the r.h.s. of Eq. (14.101) would be zero, and then the only solution consistent with 
homogeneity would be (e^) = 0, as expected in QLFT The source term explicitly breaks 
down conformal symmetry and thus introduces a mass scale in the problem, to be physically 
identified with the average value of the saturation momentum. To better appreciate the role of 
the source term in that sense, imagine starting with two different scales, say, M 2 (instead of Ql) 
in front of the Liouville potential in Eq. (14.31) and Ql in front of the source term. After shifting 
the field as = + O , with <t0 o = \n(Ql/M 2 ), the scale M 2 gets replaced by Ql in the action 
for 0. This 'scale transmutation' is generic in relation with LFT: since the Liouville potential 
is by itself scale-invariant, the associated mass parameter (M 2 in the above discussion) adjusts 
itself to the scale introduced by the symmetry-violating term, if any. 

Eq. (|4.1 II) is remarkable also in a different respect: it shows that, in the present theory, the 
one-point function (e' 7 *) is finite (and equal to one) without any ultraviolet renormalization. 
This is remarkable since, a priori, one would expect this quantity to be dominated by the fluc- 
tuations with the highest momenta and thus be afflicted with ultraviolet divergences. This was 



17 



already the case in the free field theory, cf. Eq. (13.61 ), and this is also the general situation in 
the known two-dimensional field theories, which are superrenormalizable, but not strictly fi- 
nite. The finite result in Eq. (14.1 II ) anticipates a more general property of the effective theory 
(14.91 ), that we shall verify in Sect. [5] via explicit perturbative calculations up to two-loop order : 
Namely, this theory is ultraviolet finite, in the following sense: all the n-point functions (with 
n > 2) of the Liouville field <fi, as well as all the n-point correlation functions (with n > 1) of 
the vertex operator V = e CT *, come out truly finite when computed in perturbation theory and 
thus do not require ultraviolet renormalization. (More general operators, however, like e * with 
a / a, can still meet with UV divergences, which then can be renormalized in the standard 
way.) As we shall discover in Sect. [5l this UV-finiteness comes out as a result of order-by- 
order cancellations between divergent, 'tadpole', diagrams, which are abundantly produced by 
the perturbative expansion, but which precisely cancel with each other, due to the special sym- 
metry factors of the interaction vertices. An important consequence of such cancellations is that 
the vertex operator has no anomalous dimension. 

But whereas the perturbation theory for the effective theory is meaningful and rather straight- 
forward (see Sect. [5]), it is on the other hand difficult to derive firm results about the non- 
perturbative behaviour of the theory in the interesting regime at high energy, or strong coupling, 
(j > 1 (with the noticeable exception of Eq. ( 14.1 II) ). Unlike in the standard Liouville theory, 
here one cannot rely anymore on conformal symmetry to deduce, or at least constrain, the gen- 
eral form of the correlations. In what follows we shall attempt to deduce some general properties 
of the theory from a qualitative analysis of its action (14.91) . 

Since this action involves a physical mass scale Ql, which is moreover the curvature of the 
potential at its minimum, it is quite clear that the correlation functions in this theory should 
die out exponentially over sufficiently large distances R ^> I/Qq. However, the potential in 
Eq. ( 14.91) is not a standard mass term for 0, and this difference has interesting consequences: 

First, the minimum of the potential, which roughly speaking indicates the most probable 
value for the saturation momentum, occurs at a harder scale in Eq. (14.91) than it was the case 
in the free field theory (13.21) , or in the Gaussian approximation (12.11) . Indeed, in Eq. (14.91 ) this 
minimum corresponds to Q 2 S = Ql = Ql(e a ^), whereas in Eqs. (13.21 ) and (12.11 ) it rather corre- 
sponds to Q 2 S = Q 2 = Ql e ffW (recall Eqs. (Q-Q)). Since (e^) is always larger than e aW , 
it is clear that, in the theory with action (14.91 ), the fluctuations in Q 2 are pushed towards harder 
scales, as anticipated. 

Second, unlike the quadratic potential in Eqs. (13.21 ) or (12.11) , the one in Eq. (14.91) is asymmetric 
under </> — > — </>, and this asymmetry is very pronounced for the relatively strong fluctuations 
with cr0 > 1: the potential favors large negative fluctuations as opposed to large positive ones. 
One may think that the exponential piece of this potential totally forbids the fluctuations with 
cr0 ^ 1 (or Q 2 ^> Ql), but this is actually not so: arbitrarily hard fluctuations with Q 2 S ^> Ql are 
still allowed, because they have tiny sizes and thus give small contributions to the action. Such 
a propensity towards small-size fluctuations is of course natural in any field theory, but this is 
not properly taken into account by the coarse-graining approximation (12.11) , which ignores any 
information about the sizes of the spots. 

To be (slightly) more quantitative, let us estimate the contribution of a fluctuation with size R 
to the action. The typical gradients for this configuration are V ~ 1/R, hence the kinetic term 
contributes /r(V 1 0) 2 ~ R 2 (V l (f)) 2 ~ 4> 2 . After similarly estimating the potential term, one has 
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S(R) ~ 2 + R 2 ^ (e CT ^ 



C70). 



(4.12) 



Roughly speaking, the allowed fluctuations are those for which S(R) ^ 1. 

Consider first the situation in pure Liouville theory, i.e., without the source term; then, what- 
ever the value of R is, the exponential potential allows for all the fluctuations with ^ ma x, 
with a0 max = ln(cx 2 / R 2 Ql). Note that cr0 max can be arbitrarily large, as anticipated, provided 
R is correspondingly small. The correlation functions of the vertex operator are dominated by 
the fluctuations with maximal field strength, which implies, e.g., 



in qualitative agreement (in so far as the independence is concerned) with the correct result, 
Eq. (|4.7I) . Note also that the saturation momentum for the relevant fluctuations is correlated to 
their size, Q^e "^* ^ 1/R 2 ,'m agreement with the uncertainty principle. 

We now turn to the full potential, including the linear source term. Then, we need to distin- 
guish between two kinds of fluctuations — small-size and large-size — , according to the value 
of the ratio R 2 Ql/a 2 : 

(i) When R 2 Ql <C a 2 , we expect the same situation as in pure Liouville theor\f"1. Indeed, 
in this case, the potential in Eq. ( 14.121) authorizes fluctuations within the relatively wide range 
0min ~ 4> ~ 0max> with max as above and cr0 m i n = —o 2 jR 2 Q\. The correlations of e "* are 
controlled by the fields towards the upper limit, where the source term is negligible; that is, 
they are determined by the Liouville piece of the action, and thus are power-like, with the same 
powers as in LFT 

(ii) On the other hand, for R 2 Ql ^ a 2 , the potential allows only weak-amplitude fluctua- 
tions, such that a\<p\ ^ a 2 jR 2 Q\ < 1. (By itself, the exponential piece of the potential would 
also allow for larger negative values, but these are suppressed by the source term.) Within this 
range, the potential reduces to a quadratic mass term with mass Q\. Hence, the large-size fluc- 
tuations die out exponentially over a typical distance 1/Qo- 

To summarize, the (admittedly crude) estimates above suggest that, in the effective theory 
with action (14.91 ), the correlations of the vertex operator e CT * have a power-law behaviour, with 
Liouville-like exponents, over short distances R < l/Qo, but they decay exponentially over 
larger distances R^> 1/Qq. Note that, in order to conclude in favor of Liouville-like exponents 
at short distances, it was essential that the vertex operator has no anomalous dimension, as we 
shall check via perturbative calculations in Sect. [5J 

4.3 A more general source term 

A physical hadron is never homogeneous, and the original inhomogeneity at low energy gets 
transmitted to, and it is modified by, the high energy evolution. Assume, e.g., that one starts 
with a large nucleus at Y = 0, in which the gluon density is large and quasi-homogeneous 
inside a large disk of radius R, but it rapidly drops out to zero at impact parameters larger 
than R (so like in the McLerran-Venugopalan model [36]). Whatever was the initial law for 

11 We implicitly assume here that the coupling is sufficiently strong: a <; 1. 



p4>{x) a-cf>(y) 




where 



x — y\ = R, 



(4.13) 
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this fall-off at large distances, it will get replaced by a power law after evolving the system 
to sufficiently high energies according to the perturbative evolution equations in QCD. The 
BK equation (or, equivalently, the BFKL equation supplemented with a saturation boundary 
condition) predicts that, for the saturation momentum, this power must be 4: Q 2 (x) ~ l/x 4 for 
x 3> R (see, e.g., [45,46]). It is therefore interesting to consider generalizations of the effective 
theory introduced in the previous subsection which allow for a inhomogeneous source term. 
The action then becomes S = Sl + f x j(x)(fi(x), with j(x) describing the strength and the 
impact parameter dependence of the average saturation momentum. 

The functional form of j(x) is in principle fixed by the underlying evolution equations (to- 
gether with the initial conditions at low energy) and represents a 'free parameter' from the 
perspective of the effective theory. Here, we shall parameterize this in the form 



Q 



2 

j(x) = --e* w , (4.14) 



a 



where Qq(Y) determines the value of the average saturation momentum at the center of the 
hadron and at rapidity Y, whereas the function 4>o(x), with </>o(0) = 0, describes the profile of 
(Q 2 (x)) in impact-parameter space. As previously mentioned, a physically motivated choice is 



( X ) 



[1 + x 2 /R 



212 



(4.15) 



with x = \x\ and R the hadron radius at Y = 0. The current in Eq. (|4.14| )-( |4.15l ) is clearly only 
an approximation, as it implicitly assumes that the evolutions in Y and in x decouple from each 
other; yet, this approximation captures the salient features of this evolution, namely the fact that 
the (average) saturation momentum rises rapidly with Y at any x and it develops a 1/ x 4 -tail at 
distances x ^> R. We shall furthermore assume that Q 2 I (Y)R 2 ^> 1 at any Y, as appropriate for 
a sufficiently large nucleus which is relatively dense already at Y = and has a size R fixed by 
the non-perturbative, soft, physics. 

The first Dyson equation corresponding to this current, that is (cf. Eq. (14.101 )) 

_ V 2 (0(a;)) + ^ (e^ x) ) = 9l e ^) ? (4 . 16) 



cr \ ' a 

cannot be exactly solved in general, but it clearly implies {Q 2 (x)) m Q^Q a ^ x \ (Indeed, the 
scale of the inhomogeneity being fixed by R, one has |V 2 (0) | ~ R 2 , which is much smaller 
than Ql/o at any Y .) 

It is furthermore quite clear that the short-distance behaviour of the theory (on distance scales 
r <C R) cannot be changed by the soft inhomogeneity visible in Eq. (14.151) . In particular, the 
relevant correlation functions (those of the Liouville field and of the vertex operator e a ^ x ^) are 
still ultraviolet finite. Also, these correlations preserve the same behaviour on short (r < 1/ Q ) 
and intermediate (1/Q <^ r <^ R) distances as in the homogeneous case. On the other hand, 
on very large distances r ^> R, the correlation should show a slower decay, although still 
exponential, because the effective mass for this decay, namely (Q 2 (x)), becomes smaller and 
smaller when increasing the distance from the origin. 

The role of (Ql(x)) as an effective mass becomes manifest on the expansion of the action 
S around its saddle point at a — >• 0, as appropriate for the purposes of perturbation theory. 
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Specifically, the saddle point condition is obtained by removing the brackets in Eq. (14.161) and 
determines the classical solution (j) c \{x) : 

- V 2 <t>Jx) + 5° e°**W = 51 . (4.17) 

a a 

For instance, for the particular current (|4.15l) . the above equation can be easily solved to give 

By separating = c i + 50 and expanding around c i, one finds the action which governs the 
dynamics of the fluctuation field 5(j): 



S[<i>) = S[M + J d 2 x |^(V^0) 2 + (e rf * - o5<P - 

= sm + / d 2 * {\(y*5<pf + l -Ql{x)5<? + q 2 ^) £ ^p^l > ( 4 - 19 ) 



where Q^cc) = Qo e<T * s me value of (Q 2 s (x)) in the saddle point approximation, and plays 
the role of a point-dependent mass for 50, as anticipated. The expansion in the second line of 
Eq. (14.191) will be used in the perturbative calculations to be presented in the next section (for 
the homogeneous case O = 0, for simplicity). 



5 Ultraviolet finiteness of the correlation functions 



In this section, we shall demonstrate via explicit calculations that some interesting classes of 
correlations generated by the effective theory are ultraviolet finite in perturbation theory. Our 
calculations will be performed only to finite orders (namely, up to two-loop order), and thus 
they cannot be seen as a complete, and even less rigorous, proof in that sense. Yet, as we shall 
see, they reveal a very non-trivial pattern of tadpole cancellations, which is very unlikely to be 
accidental, but most probably is representative for the way how such cancellations proceed to 
all orders in the perturbative expansion. 

The correlations that we shall show to be finite are the n-point connected correlation func- 
tions of (j) with n > 2 and all the n-point correlation functions (with n > 1) of the vertex 
operator V = e CT< ^. The mean field (<p), on the other hand, appears not to be finite, as its pertur- 
bative expansion starts with a divergent contribution of O(a) (while the respective corrections 
of higher orders are still found be finite). But this unique divergence plays an essential role in 
that it cancels, via its iterations, similar divergences which appear in the expansion of the vertex 
operator and in the disconnected pieces of the n-point functions of 0. 



5.1 The n-point functions of the Liouville field 



We shall first consider the n-point functions of the Liouville field 0, which are interesting not 
only by themselves, but also as ingredients in the respective calculations for the vertex operator, 
to be presented in the next subsection. 
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2 x y 2x 

Fig. 3. The four diagrams contributing to the 2-point function (4>(x)4>(y)) to 0(a 2 ), in the order in 
which they are listed in Eq. (15.6I ). 



As a representative example, we shall consider the perturbative expansion for 2-point func- 
tion (<t>(x)(f>(y)) up to 0(a A ), meaning two-loop order for the self-energy. Our presentation 
will focus on the systematics of tadpole cancellations, which is our main concern here. 

The expectation value is defined by the following path integral: 

(00)0(2/)} = i J D\4>] 00)0(2/) e-*M , (5.1) 

where S is the action in Eq. (14.91) and Z denotes the partition function: Z = f D [0] exp { — 5 [0] } . 
For the purposes of the perturbation theory, we separate the quadratic and the interaction parts 
of S in the standard way: S = S + S int , where 

S o [0]ee y^x^f + Q 2 ^ 2 }, (5.2) 
provides the free propagator D (that is, Eq. (13.31) with Q 2 — > Qfy, and 

=$Jd 2 x [e°* - 1 - <70 - l(a0) 2 } = J d 2 xJ2 , (5.3) 

provides the interaction vertices. The special symmetry factors associated with these vertices, 
as generated by the expansion of the exponential, lie at the heart of the ultraviolet finiteness to 
be demonstrated below: at any given order in perturbation theory, tadpoles produced by various 
vertices fc cancel with each other, because of these special symmetry factors. 

The following identity, involving the free propagator, will be also useful in what follows : 

r r d 2 k e ik ( x - z ) i 
D {x -z)^JizJ wj - 2WT ^= m . (5.4) 

By expanding 5 int in Eq. (1511 to 0(a 2 ), one finds (<f>(x)<f>(y)) « D(x, y) + (0(cc)0(t/)) (2) , 
with (an upper index on an expectation value indicates the order in a to which the respective 
expectation value is to be evaluated) 



(0(aO0(l/)> (2) = J (0(*)0(?/)0 4 (*: 







a 2 Qf 



+ w J Z1 Z2 \tt x MvWM* M\ ' (5 - 5) 

where (• • -) ee (• • -)'- '' is a simpler notation for the expectation value computed with the 
quadratic action So- After performing the above contractions, one obtains the following four 
terms (recall that the vacuum diagrams are eliminated by the denominator Z in Eq. (15.11) ) 
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Fig. 4. The double-tadpole diagrams contributing to the self-energy at two-loop order. 



(<f>(x)<f>(y)) 



= -^fi?(0) J z D(x-z)D(z-y) 
+ ^ £>(0) / -D(a? - *0Z>(*i - y)£>(«i - z 2 ) 




+ [ D(x- z x )D\z x - z 2 )D(z 2 - y) 



(5.6) 



corresponding to the four diagrams exhibited in Fig. [3l The first two diagrams involve the di- 
vergent tadpole -D(O), but the corresponding symmetry factors are such that these diagrams 
precisely cancel with each other. (Note that the integral over z 2 in the second term can be per- 
formed using Eq. (15.41) .) The third diagram yields a finite contribution in d — 2. Finally, the 
last term in Eq. (15.61) . which involves a tadpole squared, is recognized as the disconnected piece 
((f) (x)) (4>(y)) of the 2-point function (cf>(x)<f)(y)) (to this order). One has indeed: 



To summarize, the self-energy to one-loop order is given by the three connected diagrams 
in Fig. |3] (after amputating the external lines); the first two diagrams are divergent but cancel 
with each other, while the third one is finite and represents the net one-loop contribution. Note 
that the cancellation of the UV divergences has occurred between tadpoles generated by two 
different interaction vertices: (f> 3 and 4 . 

When moving to 0(a A ), the pattern of tadpole cancellations becomes significantly more 
complicated. The one-loop subdivergences (as associated with self-energy insertions on the in- 
ternal propagators) cancel with each other as explained above, but there are six self-energy dia- 
grams which feature genuinely two-loop divergences, that is, which are proportional to D 2 (0). 
These diagrams are displayed in Fig. HI together with the respective symmetry factors. Each 
such a diagram equals QlD 2 (0) times the corresponding symmetry factor. Thus, as one can 
read off Fig. HI the two-loop divergences exactly compensate with each other. The net contri- 
bution to the self-energy to this order is given by the diagrams exhibited in Fig. [51 all of them 




(5.7) 
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Fig. 5. The UV-finite two-loop diagrams which give the net result for the self-energy at 0(a 4 ). 

O $-© 4^0 



i^XD -I— CKO 



4 ^ O 8 "^^^ 

Fig. 6. All the diagrams contributing to the mean field (<j)(x)} up to 0(a 3 ). The first line illustrates the 
three diagrams which survive in the net result, cf. Eq. (15.8I) . The diagrams in the second line cancel with 
each other, and similarly for those in the third line. 

being finite. (We do not display this final result, since not particularly illuminating.) Note that, 
at this level, the compensating tadpoles are produced by vertices <p k with k ranging from 3 to 
6. Clearly, the fact that the symmetry factor for the vertex (j) k has the specific value l/k\, as 
generated by the expansion in Eq. (I5.3I) . was crucial for the success of these cancellations. 

It is furthermore interesting to notice that the two-loop contribution to the average field 
(cf)(x)}, which is of 0(a 3 ), is UV-finite as well — all the contributing tadpoles mutually cancel, 
as shown in Fig. |6] We display here the net result to this order, since this will be useful later on 
(see Fig.[6]for the respective diagrams) : 

(ftx)) = -°- D(0) + ^ / D(x- Zl )D\ Zl - z 2 ) 

/ D(x - Zl )D(zi - z 2 )D( Zl - z 3 )D 2 (z 2 - z 3 ) + 0(a 5 ) . (5.8) 

The apparent a?-dependence of the 0(cr 3 )-terms is only illusory: by homogeneity, the results 
of the integrations over z 2 in the first such a term, respectively over z 2 and z 3 in the second 
one, are independent of Z\\ hence, the integral over Z\ can be explicitly done with the help of 
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Eq. (15.41 ), and then the ^-dependence disappears indeed. 

It becomes more and more tedious to extend such explicit calculations to higher orders in 
a, and we shall not attempt to do so. However, we are confident that the pattern of tadpole 
cancellations demonstrated by the explicit examples above, as well by those those to follow in 
the next subsection, is non-trivial enough not to be accidental, but rather it is representative 
for similar cancellations taking place to all orders. Based on that, we conjecture that all the 
n-point functions of <fi with n > 2 come out finite to all orders in perturbation theory in this 
effective theory. Moreover, it seems that even for the 1-point function ((f)), the perturbative 
corrections are finite beyond one-loop order: the only UV divergence in this theory seems to be 
the lowest order (one-loop) contribution to the mean field, cf. Eq. (15.71) . From Sect.[2l we recall 
that (<p) enters the definition of the 'average saturation momentum' in the sense of Eq. (12.21) : 
Q1 = Qle a( ~^. Hence, by choosing a physical value for the latter, one could in principle fix the 
value of the ultraviolet cutoff in the effective theory. However, this appears to be superfluous 
within the present context, where Q 2 S plays no special role and thus needs not be introduced. 
Rather, it is natural to define the saturation momentum in terms of the vertex operator, whose 
correlation functions will be discussed in the next subsection. 



5.2 The n-point functions of the vertex operator 

Consider now the n-point functions of the vertex operator V(x) = e a ^ x ^ — hence, of the 
saturation momentum Q 2 s (x) = QqV(x) — , which are defined as 

(V( Xl )V(x 2 ) ■ ■ ■ V(x n )) =^jD[<f>] V(x x )V(x 2 ) ■ ■ ■ V(x n ) e- s[tt>] . (5.9) 

It is straightforward to construct Dyson equations relating these quantities to the correlations 
functions {<f>(xi)(/){x2) • • • 4>{x n )) of the Liouville field, and thus deduce the ultraviolet finite- 
ness of the former from the corresponding property of the latter, as discussed in the previ- 
ous subsection. We have already seen this on the example of the 1-point function for which 
Eq. (14.101) implies (V(x)) = 1. For the 2-point function, one similarly finds 

V* VjM*My)> + ^ (V(x)V(y) - 1) = (Ql + Vl)5^{x - y) . (5.10) 

But since Dyson equations are often formal (precisely because of ultraviolet divergences), it is 
still instructive to verify the ultraviolet-finiteness via some explicit calculations in perturbation 
theory. Below, we shall do that up to two-loop order for the 1-point and the 2-point functions 
of the vertex operator. As we shall see, the pattern of tadpole cancellations is now even richer, 
because it extends to the additional tadpoles generated by the expansion of the vertex operators. 

(a) One-loop order 

For the 1-point function (V(x)), the corresponding calculation is straightforward (as before, 
an upper script on an expectation value indicates the order in a to which that expectation value 
must be evaluated) : 
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2 

(^)) = l + ^(^)) (1) + y(0 2 (^)) (O) + 



2 

2 2 

l- a -D{Q) + a -D{0)+ ... = l + 0(a 4 ), (5.11) 



where we have also used Eq. (15.71) . But although very simple, the above calculation illustrates 
an important point which is generic: when computing the n-point functions of V, tadpoles 
generated by the interaction vertices from the action — in Eq. (15. 1 lb . the trilinear vertex which 
is implicit inside ((f)}^, cf. Fig. [6] — cancel against other tadpoles generated by the vertices 
produced when expanding V — in Eq. (|5.1 II) . the quadratic vertex explicit in the second term 
there (see also Fig. [8]). 

For the 2-point function, one needs to work a little harder. Note first that, to 0(a 2 ), (V x V y ) 
is automatically finite in the present theory — at variance to what happens in the free theory, cf. 
Eq. (13.51) — , because the only divergences that could appear to that order are those associated 
with the renormalization of the individual vertex operators, but these cancel out according to 
Eq. (15.111) . This is a general property: since we expect (V) = 1 to all orders, in the calculation 
of (V x V y ) we can discard all the diagrams contributing to the disconnected piece (V x )(V y ), 
anticipating that they sum up to one. Then, to 0(a 4 ) — corresponding to one-loop order for 
the connected piece — one finds (with V x = V{x), etc) 



iy.v v ) = l + <7 2 (d xv + (<t>^ v f>) + "- (fy v + 4>*<t>l) m 



where one should keep only the 'connected' pieces of the various correlation functions appear- 
ing in the r.h.s., that is, the contributions in which the two external points x and y are connected 
with each other. E.g., in 

£#) = 2D\x-y) + D\U\ (5.13) 

one must discard the second, divergent, piece, D 2 (0) = (<^) (<^) , since this apart of (V x )(V y ). 
With this rule, the other correlations appearing in Eq. (15.121) are evaluated as 

%4> y ) o = 3D(0)D xy , (5.14) 

and, respectively, 



-aD(0)D xy -aQ 2 f D 2 xz D zy , (5.15) 

J z 



where we have also used Eq. (15.41) . The quantities in Eqs. (15.141) and (15.151) represent one-loop 
vertex corrections at x (there are, of course, similar corrections at y; see Fig. [7]) and involve 
divergent tadpoles. However, when inserted into Eq. (15.121) . the symmetry factors are such that 
these tadpoles cancel with each other. Once again, one of these tadpoles has been generated 
via contractions inside the vertex operator V x (the one in Eq. (15.141 )) and the other one, via 
contractions inside a vertex from the interaction piece of the action (that in Eq. (15.151) ). 
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Fig. 7. All the diagrams contributing to the 2-point function (V x V y ) up to 0(a 4 ). A vertex with a blob 
represents a factor of a coming from the expansion of the exponential operators. The diagrams in the 
first line are finite and contribute to the net result, Eq. (15 - 16b - The other diagrams are divergent but cancel 
with each other. 



Finally, Eq. (15.121 ) involves the connected part of the 2-point function (0 £C y / 2 ' ) , which has 
been previously shown to be UV finite (this is given by the third term in the r.h.s. of Eq. (15.61) ). 
All the diagrams contributing to the r.h.s. of Eq. (15.121) — including the divergent ones which 
mutually cancel — are displayed in Fig. [7] By adding the previous results, one finds the follow- 
ing, finite, result for (V x V y ) to 0(a 4 ) : 



(V x V y ) = 1 + a 2 D 



xy 



a 



2 D 



xy 



e Ql 



f D(x-z 1 )D 2 (z 1 -z 2 )D(z 1 

J Zl .zo 



2 J Zl,Z2 

4/02 



y) 



o 4 Q 



D 2 D + D D 2 

^xz^zy I ^XZ^zy 



0{a% 



(5.16) 



All the terms in the r.h.s. of Eq. (15.161) except for the last one can be recognized as the expansion 
of exp{a 2 G xy }, with G xy = ((f)(x)(j)(y)) the full propagator, to the order of interest. However, 
the presence of the last term, which represents vertex corrections, shows that such a simple 
exponentiation does not hold in the present theory, in contrast to what happens in the free theory 
(recall Eq. (1X81) ). 

(b) Two-loop order 

Although considerably more involved (especially for the 2-point function), the two-loop 
calculations are also more interesting, in that they show a richer pattern of tadpole cancellations. 
We start with the 1-point function, for which 'two loops' means 0(a 4 ). One has 

(V( X )f> = a + y (^(x))' 2 ' + £ (*>(«)) W + £ , (5.17) 

where the first two terms are already known, cf. Eq. (15.81) and, respectively, Eq. (|5.6I) — within 
the latter, one has to take x = y, thus yielding 
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Fig. 8. All the diagrams contributing to the 1-point function (V^) up to two-loop order. Diagrams within 
a same group cancel with each other. 



— , while the remaining two are easily computed as 



J &(x-z) - ^D\0), (5.19) 



and, respectively, 



£<#*(«)) m = £fl»(0). (5.20) 

The Feynman graphs associated with these various contributions are illustrated in Fig.[H where 
the respective cancellations are also indicated. By adding all the contributions, it is clear that 

(V(x)) = 1 + 0(a 6 ) . (5.21) 

Consider now the two-loop contributions to (V x V y ), which count to 0(a 6 ). One of these 
contributions is a 2 ((f>(x)(j)(y)y\ which has been already argued to be UV-finite in Sect. 15.11 
The other ones involve various vertex corrections, which separately develop a large number of 
(single or double) tadpoles, but which add together to a finite result. The respective cancella- 
tions show a pattern which is similar to, but richer than, the one already observed for the 2-loop 
self-energy in Sect. 15. II In Fig.|9l we display the diagrams involving double tadpoles (together 
with the corresponding symmetry factors), grouped in such a way to illustrate the various can- 
cellations. Note that there are three types of such diagrams, each of them involving three graphs 
whose symmetry factors are such that they mutually cancel. We shall not display the diagrams 
involving a single tadpole which cancel among each other, nor the finite ones which yield the 
net result at this order, since these diagrams are quite numerous and the result is not particularly 
illuminating. Suffices to say that the respective UV-finite diagrams exhibit the same two-loop 
topologies as the self-energy diagrams shown in Fig. [5l but which now appear also in the form 
of vertex corrections. 
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Fig. 9. Divergent two-loop diagrams contributing to (V x V y ) which involve double tadpoles. The net 
result is zero on each of the three lines. 



6 Conclusions 

In this paper, we have proposed an effective scalar field theory for the distribution of the 
saturation momentum in the two-dimensional impact-parameter space, in QCD at high energy 
and large iV c . This theory has been constructed as the minimal field-theoretical generalization 
of the 6-independent distribution in Refs. [1,2] which is local (in both b and the rapidity Y) 
and which is consistent with the uncertainty principle and the conformal symmetry of the high- 
energy evolution in perturbative QCD to leading order. 

The effective action consists in two pieces: one which is universal (in the sense that it is 
uniquely fixed by our general assumptions) and is recognized as the conformally-invariant Li- 
ouville action, and the other which explicitly breaks down conformal symmetry, thus introduc- 
ing a physical scale for the saturation momentum, and which is less universal — in the sense 
that, first, there are many technical options for breaking down conformal symmetry (here, we 
have chosen a particularly simple one: a source term linear in 0) and, second, the structure 
of this piece also depends upon the details of the average gluon distribution that one needs 
to reproduce. In its minimal version, the theory is characterized by only two parameters, which 
are energy-dependent and homogeneous: the expectation value Ql(Y) = (Q 2 S ) of the saturation 
momentum and the coupling constant a(Y) which characterizes the disorder introduced by fluc- 
tuations, and corresponds to the 'front dispersion' in the previous literature [1,2]. But additional 
parameters can be added, if needed (e.g., in order to describe a inhomogeneous distribution), 
via the symmetry-breaking term. 

The description offered by this effective theory is certainly crude and oversimplified as com- 
pared to the original QCD problem. It has some obvious shortcomings — e.g., it cannot accom- 
modate the long-range correlations in impact-parameter space, which are non-local in rapidity 
— and there might be some other, less obvious, ones, which are however difficult to recognize 
in the absence of explicit QCD calculations (like solutions to the Pomeron loop equations). On 
the other hand, we expect this theory to be correct at least at short distances, of the order of 
the average saturation length 1/Q (Y) or smaller, and this is important, since this is precisely 
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the range where the predictions of the theory are more interesting, and also universal, since 
determined by the Liouville piece of the action alone. 

Specifically, there are two interesting predictions: (i) The short-range correlations are power- 
like, with the powers determined by the natural dimension of the squared saturation momentum 
Q 2 , which is 2. This result is a direct consequence of the conformal symmetry, and it is in- 
deed natural to find manifestations of this symmetry on distance scales which are shorter than 
any symmetry-breaking length scale in the problem (here, 1/Qo). (ii) The fluctuations in the 
saturation momentum are pushed to much harder scales than expected from previous analy- 
ses: the minimum of the potential, corresponding to the most probable value for Q 2 , occurs for 
Q 2 = (Q 2 S ) = <3o(e CT *), rather than for Q 2 = Q 2 = Q 2 e aW , as it was the case for the coarse- 
grained distribution in Eq. (12.11) . This is a reflection of the uncertainty principle, together with 
the multivalence of the saturation momentum: Q 2 (x), which is an operator in the present theory, 
is also the measure of the size of the fluctuation at x, and hence the natural value for its gradient. 
Whereas in the context of the coarse-grained approximation (12.11) . the value (Q 2 ) corresponds 
to a very rare fluctuation, in the tail of the distribution, this is not so anymore in the context 
of the effective field theory, where fluctuations with arbitrarily high values for Q 2 are allowed, 
since they have tiny sizes and thus give a small contribution to the action. 

It would be very interesting to understand the physical consequences of these new results. 
To that aim, the formalism needs to be further developed, to allow for the calculation of observ- 
ables. A possible direction in that sense will be sketched in the Appendix, where we show — 
within the simple context of the McLerran-Venugopalan model — how to couple the present 
field theory describing the statistics of the saturation momentum to the CGC formalism, where 
the small-a; gluons are represented by a color charge density in terms of which observables can 
be constructed in a standard way [40] . 

For our particular choice for the symmetry-breaking source term (an operator linear in 0), 
we found that the theory has a remarkable property: the correlation functions of interest (in 
particular, those of the vertex operator which enters the definition of the saturation momentum) 
are ultraviolet finite when computed in perturbation theory, and thus do not call for renormal- 
ization. Our proof in that sense is not complete (it is limited to two-loop calculations), so it 
would be interesting to extend that to an all-order proof, and also to put this result in a more 
general perspective — e.g., to clarify the general conditions under which a theory of this kind 
(a deformation of the Liouville field theory) has UV-finite correlations. 

It should be relatively straightforward to perform lattice simulations for the effective the- 
ory that we propose, and thus better understand its properties in the strong coupling/high en- 
ergy regime. The theory is finite and has no special symmetry, so there is no need for special 
care in constructing its lattice regularization. Also, lattice artifacts should rapidly vanish in the 
continuum limit. Moreover, it might be possible to analytically investigate this theory, via non- 
perturbative methods: formally, this theory is just Liouville field theory perturbed by a relatively 
simple operator — a source term linear in 0. 
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A The random McLerran-Venugopalan model 

Given the effective theory for 0, which describes the fluctuations in the saturation momen- 
tum, one can wonder how to compute the effects of these fluctuations on physical observables, 
so like the scattering amplitude of a projectile dipole. To that aim, one can resort on the CGC 
formalism, in which the observables are expressed in terms of the color charge density which 
represents the source of the gluon distribution. Using the Liouville field 0, it is possible to cou- 
ple this color charge density to the fluctuating saturation momentum in a conformally-invariant 
way. As a simple example, let us describe a 'random' extension of the McLerran-Venugopalan 
model [36], which is conformally-invariant. Recall that, in this model, the color charge density 
in a large nucleus, p a (x), is a Gaussian random field variable with zero expectation value and 
local 2-point function proportional to the saturation momentum: 

^=^|-/ ft w)' (A -" 

(J\f is a normalization constant such that fD[p] Wmv[p) = 1.) hi the original formulation of 
the model, Q 2 s (x) is a given (generally, point-dependent) quantity, which fixes the strength of 
the color charge density squared, and hence of the gluon distribution. In the present context, 
this is interpreted as Q 2 S (x) = Ql e a ^ x \ with (f)(x) a random field distributed according to the 
effective theory. This leads to a statistical field theory for the coupled fields and p a , with the 
following partition function 

Z = fD[4>) jD[ Pa ] expj-/ d 2 x Q(W) 2 + ^e^ + §|K^ } , (A.2) 

which is conformal invariant, as anticipated. To see this, note that under the scale transformation 
in Eq. (14.51) . the field p a (x), which has mass dimension two, transforms as follows: 

x -> x = \x p a (x) -> p' a {x') = —p a (x), (A3) 

so that the term coupling p to in the action in Eq. (|A.2I) is indeed invariant, so like the other 
terms there. Of course, this symmetry is ultimately broken by the source terms (for either p 
or 0), which need to be added to the action in Eq. (IA.2I) in order to generate a non-trivial 
expectation value for the saturation momentum (or for the local gluon distribution). 
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